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Renormalization group in the internal space consists of the gradual change of the coupling con- 
stants. Functional evolution equations corresponding to the change of the mass or the coupling 
constant are presented in the framework of a scalar model. The evolution in the mass which yields 
the functional generalization of the Callan-Symanzik equation for the one-particle irreducible ef- 
fective action is given in its renormalized, cutoff-independent form. The evolution of the coupling 
constant generates an evolution equation for the two-particle irreducible effective action. 

^" ' 

O ' I. INTRODUCTION 

o : 

This paper presents a modified version of the functional renormalization group (RG) method aiming at the two- 
"qj ■ particle irreducible (2PI) Green functions. The renormalization group strategy consists of the monitoring Green 
functions for the elementary field variables as the fluctuation modes are turned on gradually Q. There are two 
further ingredients which extend this method into a general purpose, non-perturbative algorithm to solve quantum 
field theories. The first is the functional generalization of the blocking where an infinitesimal fraction of the modes 
^ . are turned on or off HUSH 

This reflects the following important change of strategy. In the original version of the 
I ' renormalization group method a finite fraction of the modes are turned on or off and one relies on the perturbation 
expansion in obtaining the change in the dynamics. In order to obtain reliable results, one should work in high 
y—i • enough order of the perturbation expansion using the few bare, relevant vertices. The method of using infinitesimal 
blocking steps is based on another systematical approximation where infinitely many vertices are used in a one-loop 
• renormalization group equation. The manipulation of a large number of constants can be carried out by means of 
their generating functional. This is the reason of looking for the functional generalization of the renormalization group 
scheme. The evolution equation which looks like a one-loop equation is actually exact as long as it is solved in an 
infinite dimensional functional space. Since the general theory of functional differential equations is not well developed 
we have to rely on a truncation of the evolution equation, the only approximation of the method. It is based on the 
projection of the functional differential equation into a functional subspace. The tailoring of this subspace seems to 
be an interesting and a physically better motivated approximation scheme than the truncation of the perturbative 
series. 

In order to motivate another modification of the original strategy we recall that the order in which the modes are 
eliminated in the renormalization group method, 



o 



- 1—1 

X 



is usually given by means of a scale parameter in the external space of field theory, the space-time, such as the 
momentum or energy and the dynamics is studied as the function of the UV or IR cutoff. In this manner the short 
distance modes whose dynamics is rendered more perturbative by the kinetic energy are eliminated first and the 
more problematical nonperturbative modes are treated later when the blocked action is sufficiently renormalized, 
Pi < P2 < • • • < Pn—i 5~ Pn- An alternative to this procedure is the ordering of the modes according to their 
scale in the internal space, their amplitude. The sharp cutoff version of this procedure consists of the study of the 
A-dependence of the path integral 

Z x = f ( X # n e-* s "^<->«. (2) 

A more natural, smooth cutoff realization of this strategy is the functional generalization of the Callan-Symanzik 
equation 0, where the mass is evolved in a relativistic theory 0, , 

Z = [ D[^]e-T, S ^~^L€. (3) 



Since the mass term dominates the dynamics of free fluctuations whose wavelength is longer than the inverse mass 
the two procedures produce similar scale dependences. 
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There are in the same time essential differences between the external and the internal space realizations of the 
renormalization group. The most important one is that the physics is kept fixed along the renormalized trajectory of 
the external space method but changes together with the control parameter which suppresses the fluctuations. This 
difference arises from the different ways the UV modes are handled during the evolution. In principle both methods 
are similar, their evolution equations look formally equivalent, only the actual form of the suppression of the modes 
is different. We may introduce a straightforward division line between these methods by stating that a suppression 
mechanism corresponds to the external or internal space scheme if the corresponding evolution equation is UV finite 
or divergent. The internal space methods require the introduction of an UV regulator which is kept fixed during 
the evolution and the final products of the evolution are the bare Green functions, as opposed to the external space 
methods where the gliding cutoff renormalizes the Green functions automatically. 

In the first part of this paper we present a modification of the internal space evolution equation which eliminates this 
shortcoming and yields UV finite, renormalized Green functions. This is achieved by the functional implementation 
of a renormalization which is reminiscent of the BPHZ method [Tol ITU IT^ | . One finds as byproduct a possibility of 
a Hartree-type optimization of the unavoidable truncation of the evolution equation. These issues are studied in the 
framework of the effective action for the field variable, the generating functional of the one-particle irreducible (1PI) 
Green functions and serve as a preparation for the second part of this work, the construction of an internal space 
renormalization scheme for the computation of the 2PI Green functions. 

This step is made by generalizing the Callan-Symanzik scheme into a far more flexible procedure by allowing the 
control parameter, the mass square in the preceding case, to be replaced by any parameter of the dynamics which 
controls the physics in a differentiable manner. Such kind of an extension is needed when symmetry considerations 
require the presence of soft modes and the suppression of the fluctuations by a mass term is inconsistent. Actually, the 
mass controls the amplitude of the fluctuations but disregards their nature. The Gaussian fluctuations in the absence 
of interactions are rather peculiar, their higher moments factories according to Wick's theorem. The interactions 
introduce non-factorisable correlations. 

It is easy to identify a control parameter which turns the non-factorisable fluctuations on. In fact, it is just the 
coupling constant or a common multiplicative factor in front of all coupling constants if there are several of them. 
The evolution equation, generated by the coupling constant is highly involved |8(. In fact, the evolution equation 
contains one-loop integrals as long as the suppression term which generates the evolution is quadratic in the field. The 
complications arising from a quartic suppression term, </> 4 , the two- and three- loop integrals in the evolution equation 
reflect the complicated structure of excitations created by the operator <j> in terms of single-particle excitations. One 
can simplify the situation, at least formally, by using composite operators in which the suppression term is again 
quadratic. The evolution equation for the effective action for such composite fields will be as simple as that of the 
Callan-Symanzik scheme. In the case of the </> 4 scalar theory the composite operator one turns at is obviously 4> 2 
and one expects that the gradual turning on the interactions should be possible to realize in a simple manner by 
following the evolution of the effective action for the composite field (j> 2 . But the physics encapsulated by such an 
effective action is far more restricted than those found in the effective action for the non-local composite operator 
4>x4>y |l3lll4j . In fact, the polarization cloud around a bare particle requires the possibility of separating the two local, 
elementary field variables in 4> 2 in order to trace the space-time structure of dressed, physical excitations. Therefore 
our proposal, spelled out in the second part is to describe the structure of vacuum up to two-particle correlations by 
solving the evolution equation for the effective action for the propagator. 

This scheme is useful if the initial conditions, the effective action at the initial value of the coupling constant is 
known. One should therefore give the initial conditions for vanishing or infinitesimally weak coupling strength. The 
evolution equation will then be integrated until the physical value of the coupling constant is reached. We shall 
investigate in this work the solution arising from the non-interacting initial condition only. Such a scheme has already 
been found useful in performing gauge- free computations in QED |l5j . in constructing systematically the density 
functional ^(| , and for identifying the localization-localization transition for disordered systems 01 • 

There are formal similarities between the Schwinger-Dyson hierarchy of equations and the evolution equation. But 
the main advantage of the evolution equation approach is that it provides a coupled set of equations for the Green 
functions of the composite field which are uniquely solvable for a given initial condition in the perturbative regime. 

The organization of the paper is the following. In Sect. [H]a short description of the functional Callan-Symanzik 
approach is given, its solution by iteration and the improvement of the truncated iteration by optimization is pointed 
out. Sect. IIIII presents the derivation of the internal space RG evolution equation for the 2PI effective action for a 
family of composite fields, the bilinear expressions of the elementary field. The RG equations for the effective theory 
are presented in Sect. IIII Bl after expanding the 2PI effective action around the vacuum. Instead of their numerical 
solution an optimized iterative approach for solving the evolution equation of the 2PI effective action is discussed in 
Sect. IIII CI and the analytic solution found truncating the iteration just after the first step. The effective theory for 
the composite field is related to the theory for the elementary field in Sect. IIVI Finally, the results are summarized 
in Sect. 
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II. 1PI EFFECTIVE ACTION 



The simplest realization of the internal space renormalization group method is the control of the amplitude of 
the fluctuations by the mass term which yields a functional generalization of the Callan-Symanzik equation. As the 
simplest example we consider the Euclidean scalar model defined by the bare action 



S[<p] 



see Appendix ^ for the notations and conventions. The generating functional 



(4) 



(5) 



is given by a regulated path integral where the UV regulator is not shown explicitly. It gives rise the effective action 

6W[j] , 



-W\j]+j-<f>, 



satisfying the relations 



ST[(/>] S 2 W[j] 5 2 T[<j>] 
8(f> ' SjSj S(f>5(j> 



1. 



(6) 



(7) 



A. Evolution equation 



The dependence of W\j] and T[<j>] in the mass /i 2 , considered further on as our control parameter, is related, 

= -d^w\j\, 



(8) 



and the evolution equation can be obtained by simply bringing the derivative d^i into the functional integration in 
Eq. Q, 



d^T[4>] 



_5 2 W\j] , fSW[j] 



s 2 m 



(9) 



It is worthwhile mentioning that the harmless looking step of bringing the partial derivation with respect to the 
control parameter of the evolution within the path integral produces in general highly non-trivial results. First, it 
can be considered as the generalization of the Hellmann-Feynman theorem |l8l [l9j for time-dependent processes. 
Second, different choices for the control parameter A produce important, well-known equations: (i) When A = tp x , 
the functional derivative with respect to one of the integral variables is vanishing. The stability of the path integral 
against the shift of the integral variables produces the Schwinger-Dyson hierarchical equations (c.f. App. ID (ii) If 
A is a parameter of a transformation of ip x which preserves at least a part of the action invariant then the results are 
the Ward identities. 

Here we use the simplest realization of the RG idea to control the amplitude of the field fluctuations, namely 
the control is achieved via the mass term A = /i 2 in a momentum-independent manner 8] in order to stay as close 
as possible to the discussion of the evolution equation for the 2PI effective action presented below. There are other 
choices which allow one to control the amplitude. The most widely used suppression methods imply strong momentum 
dependence HUSH. 

The graphical interpretation of Eq. is obtained by following the Callan-Symanzik strategy. In fact, the /s 2 - 
dependence of a Feynman graph F with legs amputed arises from its internal lines, the propagators G, 



V p 



(10) 



The evolution equation is the application of this chain rule for the 1PI graphs in the framework of the skeleton 
expansion. 

Contrary to the strategy of the Callan-Symanzik equation where one expands around the massless, critical system 
our initial condition is supposed to be in the perturbative regime. Therefore we choose fly 3> A 2 where A denotes the 
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UV cutoff as initial condition because the quantum fluctuations are suppressed and the tree-level approximation to 
the effective action, 



There are two related reasons to perform further adjustments on the evolution equation ©. One is that the control 
parameter used in the internal space evolution does not regulate the UV divergences. As a result Eq. has to 
be first regulated and the resulting effective action must be renormalized in order to extract Green functions which 
correspond to cutoff independent scales. Another, more technical looking complication is that the integration of the 
evolution equation @, starting with the initial condition (|11|) . generates more and more complicated </>-dependence 
in the effective action as (x 2 is decreased. The two issues are related because once the effective action is expressed in 
terms of the renormalized quantities it is supposed to contain slower, i.e. easier to truncate dependence in the field. 

The unavoidable truncation of the evolution equation should be as harmless as possible from the point of view 
of the terms appearing in the inverse of the right hand side of the evolution equation 10 . There are two different 
directions to follow in the search of optimization. One is the appropriate choice of the dependence on the control 
parameter of the evolution in the Lagrangian. Because all we want to achieve is an interpolation between a simple, 
calculable initial condition and the true system only, different possibilities might be available to introduce the control 
parameter. This avenue has already been explored in the context of optimizing the blocking transformation in the 
traditional realization of the renormalization group method in external space [2(J, |21| . We cannot follow this direction 
here, having already chosen the mass as control parameter because of other reason. The other dimension of the 
optimization is the tailoring of the functional form of the effective action. One should recall in this respect that the 
magnitude of the variable <fi has nothing to do with the strength of the quantum fluctuations since <fi is a book-keeping 
device only in defining the generating functional for the 1PI Green functions. Nevertheless <fi is supposed to be small 
in an indirect manner for an optimized scheme because the more complicated terms of the effective action tend to be 
of higher order in <fr. More precisely, the omission of the terms O (</> 2n ) leads to ignoring the n-particle correlations. 

We shall try to slow down the ^-dependence of the effective action by the appropriate truncation of the inverse on 
the right hand side of the evolution equation. The inverse will be written as a series where terms appear in the order 
of increasing complexity in such a manner that the more complicated terms are as small as possible. The tailoring of 
the "free" inverse propagator around which the inverse is expanded is reminiscent of the Hartree approximation when 
one adds a mass term to the free Lagrangian, subtracts it from the interaction part, and optimizes the perturbation 
expansion with respect to this rearrangement of the Lagrangian. But one should notice that the same strategy 
is followed in establishing the renormalized perturbative series, too. Therefore we separate off local counterterms, 
Tct^], from the effective action, 



in such a manner that the renormalized effective action, [</>], converges as the cutoff is removed. While this 
separation has no effect on the exact solution, it will influence the solution of the truncated evolution equation. The 
condition A <C ^ clearly prevents us to send the cutoff to infinity. But an approximate cutoff independence of r^ [0] 
can still be tested for /i 2 <C A 2 <C /i 2 , for a sufficiently large but finite /j,q. 

The functional Taylor expansion involving the renormalized 1P1 Green functions T n Xl ....,x n and their counterterms, 
Tct n n,.,i„, will be used, 




(11) 



is applicable. 



B. Renormalization 



T[<t>] = r K [0]+r OT M 



(12) 




(13) 



where double indices are integrated over, c.f. Appendix The evolution equation 




(14) 
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is then written by means of expanding the inverse on the right hand side around the renormalized propagator, 
Gr = r^ 1 , and by splitting the resulting series into the sum of two pieces in d dimensions, 



d^RM = ? E (- 1 )™ Tr 



*=[#]+! 



\n=2 v 



2) 



T27J x 3 ,...,2;2 re < / ) 2:3 ' ' ' ^a:2» ' Gr 



+ \<t>'^ (15) 



and 



c^r CT [</>] = -Tr 



fs 2 r R [<f>] , <5 2 r CT [<^ 1 

2 



E (- 1 )'" Tr 



fi+i 



- 



(2ra-2)! 



2n X3,...,x 2 n < Px 3 ' ' ' < Px 2n 



Gr 



(16) 



where the first two space-time indices of r„ were suppressed. This split is made in such a manner that the renormalized 
part remains " finite" , i.e. approximately cutoff independent at energy scales p as long as p 2 <C A 2 <C Mo ■ The expansion 
of the inverse on the right hand side of Eq. around Gr, 



d^TcT^] = £(-i)-Tr < G} 



m— 

2 



Tct 2 + E (27i - 2)! X3 '---' X2n + ^ CT 2n x 3,---^2 n ) < t>x 3 ■ ■ ■ 4>x 2n j Gr 



n=2 



J2 M) m Tr G R 

=[#]+! ^ 



g (2^2)! 



2n 2J3,...,X2T 1 Y'a;3 ' ' ' < Px 2 



,Gr 



(17) 



shows that each contribution to TctI^] is UV divergent. 

The identification of the coefficients of the first few powers of (j> on both sides gives 

d^T = 0, 

^/J. 2 ^ 1 4,a;i,X2,X3,X4 — 0; 
^fi 2 ^6,xi,X2,x^,X4,x^,xe : 



2-2 3 



E Tl " [^ r 4 ^ (1) , Il(2) Gjjr 4 x x( 3 ) , a ;„ (4) G_ R r4 Xni5) ,x„ (e) GR], (18) 



and 



^r CT o = - (~l) m Tr [G K (r C T 2 G fl ) m ] = -Tr (T 2 + T CT a)" 1 , 



^sTct 2,xi,x 2 



m=0 

, oo m 

4 E(- X ) m E E Tr [ G «( r ^ 2 G fl )^ 1 G fl ,(rcT 2 G fl ) m -'''(r4 + r CT 4 )« wW> x w(a) G H ] , 

m=l 7TGS2 J = l 

. oo m — 2 m— J— 2 

sE^^EE E Tr [G fl .(r C T2G i? r--'- fe - 2 (r 4 + r CT4 )^ (1) ^ (2) G fl 

m=2 7reS 4 J = l fc=0 

x(r CT 2Gfl) : '(r4 + T C T A)x„ (3) ,x 7l{i) GR{TcT 2Gfl) fc ] 
, oo m 

+|EM) m EE 

m=l 7reS4j = l 

xTr [G R (r CT 2 G fl ) i - 1 G fl (r CT sG^)" 1- ^ + r CT Gfl] , (19) 



for the first few coefficient functions. 



6 



C. Iterative solution 



It is instructive to imagine the solution of the evolution equation as the fixed point of the iteration T^ N 



r™ + ^/. ^Tr 



''5 



2 2 

M - Mo 



/Ti 



(20) 



The speed of convergence of the iterations depends on the choice of the starting point. Let us assume in the spirit of 
the perturbation expansion that the iteration reaches at least a fixed point when the tree- level starting point, given by 
Eq. (jl 1|> , is used. The iteration in Eq. 120(1 becomes well-defined at each step only after imposing the initial condition 
pi(l together with T CT ^ [(f)] = for the /i 2 -integration, i.e. the dressing is assumed to be completely suppressed at 



Mo- 



For the sake of simplicity, we apply the local potential approximation, i.e. we restrict the field variable to be 
homogeneous, <j) x — ► $, in the effective action which will be written as 



!>[$] = V[U R „ a ($) + U CT „»(#)]. 



(21) 



The starting point is a /^-independent potential, 



U 



[o] 

R m 2 



($) = C/„ 2 ($) = + 



U, 



[o] 



\U 



(TV) 



and the iteration l|20|l becomes an integral equation, 

U CT ^ W - U CT p 2 + B CT p? a -^^ VJ R ^2 , U CT 

for the local potentials which is written in d = 4 as 



(22) 



(23) 



-2 



dp 



1 



2 



p 2 + K~^) 

# 2 E K V R^)GR^] m G 

m=3 J P 



(1 - Vr ~A$)Gr ? + \V£ p (<S>)G R ^] 2 ) G R A . 



R u 2 , 



b ct v? Q ->^[U R >2 , U CT ^3] = -- I dp 



+ (1 - Vr p (^)Gr A 2 + \Vr ^ (*)Gr ^} 2 ) G R ^ 

-\ J 2 ^ 2 J [1 - V£ , 2 ($)G fl ^ + [V R ' ^)G R ^Y] G R ^ 

ji r l d/if(-ir i - ( ^ w Y 

2 ^ 2 hh ' p 2 + u 'r^)\p 2 + u r^)) ' 



with the help of the renormalized propagator 



G R p 



p 2 + KJV 



(24) 



(25) 



and the subtracted potential 



(26) 
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The multiplicative factors h in front of the one-loop integral on the right hand sides indicate that each iteration 
resums a successive order of the loop-expansion within the local potential approximation. The result of the first 
iteration is called independent mode approximation because the quantum fluctuations are treated independently. The 
main simplification in computing in this order comes from the fact that the initial action is independent of /x 2 , 



4! 



n=0 



£$2 \» 



1 



p 2 + fl 2 J p 2 + fl 2 



M _$2 , 9 ^4 
I 



2 $ 2 + |$4 + jjU^. ^ g) _ jjU^. fa g) 



lnS±^ + ^ 2 



2^ 



1 



p \p 2 + fl 2 J p 2 + p? 
1 1 



p 2 + l4 2 \p 2 + fi 2 p 2 + fi 2 Q J 2\2 



1 ( U 2 



(p 2 + n 2 ) 2 (p 2 + fi 2 ) 2 



(27) 



where 



U ll ($;(i 2 ,g) 



. fln{p 2 + fi 2 + U 2 



stands for the one-loop effective potential of the model of mass fi and coupling constant g and 



U 1 R \^ l i 2 ,g) = U u {^,n 2 ,g)-- 



\n(j> 2 + fi ) + ^ 



9^2 1 



1 



2 p 2 + /i 2 2 \p 2 + fi 



£$2 

2^ 



(28) 



(29) 



represents its UV finite, renormalized part. When the limit fig — > oo is made then the O ($ 2 //ig) terms become 
negligible and we find the standard result, 



R v 



4! 



\n(p 2 + fi 2 ) + ^ 



9*2 1 



1 



2 p 2 + fi 2 2 V2 / (p 2 + fi 2 ) 



(30) 



up to a constant in the renormalized potential. It is easy to see that the subsequent iterations preserve the UV 
finiteness of J7r(<I>). 

It is worthwhile mentioning that the split of the iteration of Eq. I|2U|I into a renormalized and a counterterm part, 



V 1 R 



m=[f] + l 



G, 



1 

51 0„ _ 



i[N] 



(2n-2)! 2ns3,...,x 2 „^3 



+ 2>-0, ( 31 ) 



and 



- ^(-lrivc, 



i CT2 + j"2 n _ 2)! 2 ™ a5 3'---' X2 '> crsn i3,...,un™ 1 » 



m=[|]+l 



(On - 



(2n-2)! 2n ^,..,^^3 



(32) 



realizes a particular BPHZ renormalization scheme. In fact, r^[</>] is a finite functional for each N and the coun- 
tcrtcrms in rj^l[0] correspond to the full subtraction of the overall divergences at each order N which are defined by 
the expansion of the integrand of the loop integral, the trace on the right hand side of Eq. (|32|l . 



III. 2PI EFFECTIVE ACTION 



Our goal in this section is to derive the evolution equation for the 2PI effective action for the model defined by the 
action This is achieved by generalizing the standard procedure, presented briefly in section^] for the bilocal field 
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ipxfy For this end we introduce the generating functional 

together with its Legendre transform, the effective action for the two-point function, 

5W[J] 



u 2 X >V 



which is given by 



f [G] = -W[J] + ilr g f ■ J. 



The inverse Legendre transformation is based on the relation 



5Q X 



2 Jx>v 



(33) 



(34) 



(35) 



(36) 



Since the source J XtV is symmetrical, so is Q xy . Another important symmetry, translational invariance is supposed 
to be recovered in the limit J — > 0, therefore the propagator in the vacuum is taken to be translation invariant, 
(T(ip x ip y ))j = o = Q v x<y = Q v x - v - For the sake of simplicity we shall assume that the bare action is always chosen in 
such a manner that there is no condensation, (ip x )j—o — 0, as in Section ITT1 
The relation 



1 (x 1 ,X2),(yi,V2) ~ r > \ Xl,yi X 2 ,V2 "+" °Xl ,1/2 °X 2 ,Vl ) ~ 1 



cf. Eq. (| A5|) . can be used to establish the identity 

8 2 W[J] 



S\J8\J 





: I 




~6 2 r{g}~ 




sgsg 




sgsg 





(37) 



(38) 



where the inverse A~ 1! = A^ 1 : I with respect to the operator I was introduced. The matrix / acts on the two-particle 
subspace by projecting into the symmetrical subspace and plays the role of the identity operator for indistinguishable 
particles. 



A. Evolution equation 



We shall use the coupling constant g to control the RG evolution in internal space and let it go from zero to its 
physical value. The corresponding evolution of T[Q] is governed by the evolution equation 



d g f[g) = -d g W[J] = e -**M / V[v\d a S[v\e-^+^- J -* = - 



S 2 f[G] 

sgsg 



-a 



(x ,x) ,{x ,x) 



(39) 



It is worthwhile mentioning that the onset of the condensation is signaled by the vanishing of the restoring forces 
to the equilibrium position, acting on the quantum fluctuations, i.e. by g^ 1 — 0. Therefore, the exclusion of the 
condensation induces the restriction g~ x ^ 0. 

It is easy to understand the evolution equation for the functional W[J] graphically. In fact, let us assume the form 



W[J]=Y J -,WnX 1 ....,xJ Xl ---Jx n , 



(40) 



n=0 



where the repeated double-indices Xj = {xj, yj) are integrated over and the evolution equation reads as 



n=0 n=0 Jx 



W n+2 (x,x),{x,x),X 1 ,...,X n Jx l • ■ ■ Jx„ 



-ft- 



3! 



]W m+1 ( x ,x),Xi,...,X m JXi ■ ■ ■ Jx m W n+1 (x,x),Y 1 ,...X n JY 1 ■ ■ ■ Jy n 



(41) 



m.n— 
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Consider the contribution O ( J ) on both sides. On the left hand side we find d g acting on the Green function with k 
pairs of field variables which is connected as far as the cutting of two particle lines are concerned. This Green function 
is the sum of Feynman graphs and the dependence on g appears through the factors g multiplying the vertices. 
Therefore, the right hand side is supposed to be the sum of these contributions. We turn now the four internal lines 
attached to the vertex on which d g acts into external lines and consider the left-over graph. The first and the second 
terms on the right hand side stand for graphs whose left-over graph remains two-particle connected or falls into two 
disconnected components, respectively. 

Motivated by the Hartree-type optimization of the truncated functional Callan-Symanzik equation presented in 
section [H] we split off a quadratic term from the effective action, 

T[g} = T[g} + ^a XuX2 g Xl g X2 , (42) 

where a denotes a g— dependent coefficient function which is vanishing for g = 0. The corresponding evolution 
equation is 

vw-a/.[® + C M ^ +, ' ! (a i -H !R (43) 

with Lt Xi y\r UiV ) — S Xty S u>v 5 XjU . It is furthermore useful to separate the effects of interactions in the effective action 
by means of the parametrization 

r[g} = Tf ree [g} + u[g], (44) 

where the first term corresponds to the free theory, cf. Appendix iBl 

and the 'potential' U[Q] stands for the effects of the interactions. The initial condition, U[Q] — 0, for the evolution 
equation is imposed at g = 0. 



B. Expansion around the vacuum 

We return now to the evolution equation l|43[) which cannot be solved exactly. Its usefulness depends on the 
possibility of finding a good enough truncation, a projection into a restricted functional space. This latter should be, 
on the one hand, simple enough that the truncated equations can be solved either analytically or numerically and on 
the other hand, rich enough to incorporate the important effective vertices of the theory. The natural procedure is 
the functional Taylor expansion of the type Eq. (|40|l for U[G], truncated at order N, i.e. we retain the ./V-particle 
effective interactions. In order to render the Taylor series faster converging and to recover the 2PI vertex functions 
the Taylor expansion will be organized around the vacuum by splitting the propagator as g = g v + TL where g v is the 
propagator in the vacuum and writing the effective action as 

T[H] = T free [g v + H}+ 7 [W], (46) 

with 

N 1 

7 (W) = U[g v +H]=Y - 7 „ Xl ,-,x n H Xl ■■■H Xn . (47) 

The coefficient functions are symmetrical with respect to the exchange of any pairs of their variables, r y n ...,x -,...,x k ,... 

x . x, .... and are translation invariant because they correspond to the functional Taylor expansion around the 
translation invariant vacuum. This justifies the notation 

7n,Oi,«i),-,fa».8») = 8pi+qi+-+p n +q n fl%,(p 1 ,q 1 ),...,(p n ,q n ) (48) 

in momentum space. The functional Taylor expansion for the full effective action thus reads as 

oo 

r [w] = H -i r " ^ n n Xl ■■■n Xn , (49) 

^-^ n! 

n=0 
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where the first few coefficients are 

To = 

r i;0,y) = 
1 3;(xi,2/i),(x2,J/2),(a;3,2/3) — 

where 



:Tr a 



/ 2 



Tr In g v + 70 



S 2 S[0} \ 

Xi,X 2 +72.XLX2. 

X! (y„(i),s„ (2 ))^v (y„ ( 2)^,(3))^v (t/, ( 3),s„(i)) +73;(xi,yi). (^2,1/2), Osa/s)) 



2 • 3! 



2-4! 



EC-i q-i q-1 q-i 



7r£S 4 

+74 



= 6L . 



(x,y),(u,v) ^(x,y),(u,v) 5Q xy 5Q uv 



5 2 T free [g] ft _ 



see Appendix [HI The propagator in the vacuum satisfies the equation 



where the self energy is given by 

The evolution equation can now be written as 

h 



S = 2a :Sv + 271. 



(50) 



(51) 



(52) 



(53) 



9 ' T = 4! 



( 1 \ 00 1 

: ( Y\ L ~ 5 s a ) : ^ + X! fa _ fg r ";-gi,-z 2 , - ,z n Hz 2 ■ ■ -Uz^dgGw z 1 + Gv ■ a : <9 ff £ v , 



(x,a:),(x,x) 



-(n-l)! 



(54) 



by using TL as an independent variable and suppressing the first two double-indices of the coefficient functions T r , 
Here 



k x]y = x,y + 72,x,y 



(55) 



stands for the inverse of the two-particle propagator and the separated O (G 2 ) term in the effective action is treated 
by the expansion. The identification of the coefficients of different powers of TL on the two sides up to O (TL 2 ) yields 

d a % = ~J K M(X , X) + lJjK:(a: i^) ro ] (X)X)i(X)X) + G v : (~L - \d g a^j : Q v + (T t +Q V : a) : d 9 S v , 

<9 5 7i,x : 



ft . H - , 

^K Zu z 2 ^3-z 2 ,z 1 ,x + ^jTr K : (a: K : T 3 . x + r 3;X : A : a) 



L>x,zGv z - dgax.zGv z + {^2-x,z + ax,z)d g G v z, 



with 



X-, 



1 Kzi,Z2 (^3;Z 2 ,Z 3 ,X„ W ^ 3l Z 4 r3;Z4,Zi,X ir( 2 ) — ^4;Z 2 ,Z! ,X„ (1) ,X„ 

<9 ff ax 1: x 2 + T , i,x 1 ,X2,zdgG v z, 



A "x,(«,«) A "(u,tt),Y- 



(56) 



(57) 
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Iterative solution 



We follow the general strategy of Section [D] and rewrite the evolution equation iffify as integral equation, used for 
iteration, 



'^\g\ 



-a 



9 



1 



(58) 



We choose a = gaL/12 with a ^-independent constant a in the Hartree-term which guides the evolution and seek the 
effective action in the form 



The iteration of the functional fl g [Q] , 



r g [g} = a g [g] + ^^-g:L:g. 



rt»+V[Q=nW[g\ + B[nm, 



is realized by 



B\a g } = j } fdg"TrL 

' J 



S 2 n g ,[g] g'a 

sgsg + 12 



-a 



The starting point is the ^-independent effective action fig ' = Y^ ree and its first iteration, 



fl[nM] = 2Tr {[ln(g'-i 



L-i + ID-ln^rL- 1 )] : /} , 



yields 



^W^] = ^^ ee [£] + -Tr 



S l-j - l-j ■ 2 

The expansion of the logarithm on the right hand side gives 

-ga 



I : In 1 



12 



2 " n 

n=l 
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Tr [(£ : + ^-^5 ■ L : g, 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



since I : L = L. This effective action corresponds to the sum of ring diagrams where the two-particle propagators are 
joined by the vertices g. 

For the choice a = 1 the right hand side of Eq. I|58|) is O (H) and each iteration resums a successive order of the 
loop expansion to the effective action, c.f. Eq. (|53|l . This resummation explains the infinite series in r^'ft/]. When 
a = is chosen then the right hand side of Eq. I|58|l is O (g) and each iteration resums an order of the perturbation 
expansion in g only. For any other values of a an effective interaction strength is used in the partial resummation of 
the loop expansion (second terms on the right hand sides of Eqs. (|58ll and Q64[l) which is corrected order-by-order in 
the perturbation expansion (third terms on the right hand sides of Eqs. and JH3J). 



IV. RELATIONS BETWEEN FUNCTIONAL EQUATIONS 



We comment in this Section first the relation of the evolution equations with well-known, non-perturbative iterative 
methods and next the similarity of the derivation leading to the evolution equations for the 1PI and 2PI Green 
functions, Eqs. (|14f> and l|44|) . respectively. 

The evolution equations, Eqs. © and (|39f) represent exact relations for the generating functionals which, assuming 
the proper initial conditions, lead to the determination of the 1PI and 2PI Green functions, respectively. It is natural 
to raise the question about the relation of these equations to the Schwinger-Dyson functional equation which is another 
exact equation satisfied by the Green functions. Both share a common mathematical origin, t hey express the way 
the path integral for the transition amplitudes transforms during the change of certain variables |22| . The Schwinger- 
Dyson equation, overviewed briefly in Appendix [Dl is a hierarchical set of equations as the evolution equations. The 
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Schwinger-Dyson equation for the two and four point irreducible functions, Eq. (|D14|I and i|D21(l . could be converted 
into an evolution equation by acting on them by the operator d\. But the difference between the resulting equations 
and Eqs. @ and l|39(l is that the formers contain the derivative of the Green functions with respect to the control 
parameter in an implicit manner. The evolution equation scheme is organized in such a manner that the derivatives 
are all explicitly expressed in terms of the Green functions themselves. 

There are a few advantages when the evolution equation method is used. One problem with the Schwinger-Dyson 
equations is that it is not clear if its solution is unique. Their iteration may have different basin of attraction and the 
identification of the physical solution might be rather difficult. As long as the initial conditions are properly given in 
a nonphysical, but perturbative domain, the evolution equation method produces a unique solution. 

Let us now compare the evolution equations for the 1PI and 2PI functions more carefully. A remarkable property 
of the evolution equation © is that it does not lead out from the space of quadratic functionals. Due to this reason 
it is essential to keep the interaction term in the initial condition. Another way to characterize this feature is to 
notice that the structure of the evolution equation reflects the choice of the /^-dependent part of the bare action 
only. The information about the bare interaction is contained in the initial conditions to the evolution equation. The 
2PI evolution equation, Eq. Ij39(l . on the contrary reflects the O (</? 4 ) contact interactions only and is independent of 
the O ((/? 2 ) part of the bare action. In general, the evolution equation for the n-linear expression of the local field 
variable <p x contains the O (4> n+1 ) and higher order pieces of the bare action whose terms up to O (</>") appear in 
the initial conditions only. An important consequence of this structure is the UV finiteness of Eq. I|39() . In fact, the 
integrand of the formal one-loop integral on the right hand side is given in terms of the independent variable Q and as 
it stands, Eq. H39[) is a formal expression and recovers the well-known divergences for the choice Q p = O (p~ 2 ) only. 
The counterterms to the 2PI effective action depend on the choice of the independent variable Q. It was due to such 
a complicated renormalization that we left the discussion of the 2PI evolution equation on the bare level in Section 
lllll and counterterms like in Section UTI were not introduced. As of the initial conditions are concerned, they must fix 
that part of the effective action which is not represented in the evolution equation, in particular the bare interactions 
or the free part in the 1PI and 2PI cases, respectively. 

Finally, the similarity between the 1PI and 2PI evolution equations represents an opportunity to find an alternative 
to the matching strategy for the introduction of effective theories. This similarity suggests to replace the composite 
operator f x f y by a bilocal effective field ip x ,y and inquire about the corresponding effective theory given by the action 
S c s[ip]- The construction of an effective theory for local composite operators is well established by means of the 
perturbation expansion (2^ , or of non-perturbative scheme [U . Our starting point for the derivation of the effective 
theory for the bilocal field tp Xj y is the generating functional W[J] of Eq. (|33fl for the connected Green functions of 

the composite operator ip x ip v . The free bare effective action, S^[tp], is chosen in such a manner that it reproduces 
the free propagation of two identical particles controlled by ip Xl y, 

J V[ip]e~^ ^ [ ~^- vAv * + ^ m2v ^ + ^ v ' J ' v = J T>[ip]e~x s *t W+m^ :J . (65) 

The interactions can be incorporated by adding the piece g J x ip x X /A\ to the effective action. The equivalence of the 
interactive theories, expressed by the relation 

e m-J] = J D[^] c -isS/[*]-A/.<.+A^ ) (66) 



i i J '.«*7 J '.' (67) 



can be verified by acting with the differential operator 

_ hg f 

e 

on both sides of Eq. I|65|l . 

The first step of this argument, the construction of the free action S^[ip} is far more complicated than the second 
step, the inclusion of the interactions. This is not by accident, the non-interacting dynamics of equivalent particles has 
to include the symmetrization/antisymmetrization of the quantum states. This step introduces long range, topological 
correlations among the otherwise non-interacting particles. These correlations represent a far more difficult problem 
than the local interactions. This difficulty is particularly explicit in the construction of the density functional [l6| 
for non-interacting particles when the effective action has strongly non-local vertices of infinitely high order without 
natural small parameter to justify any truncation. When the non-interacting particles are fermions this procedure 
amounts to bosonisation. 

The difficulty of obtaining the free effective action is avoided when one uses the evolution equation method and 
pays the price of having an effective action for a bilocal variable Q x ,y rather than the density functional which depends 
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on the local, diagonal part Q x .x only. The complications mentioned are then hidden in the logarithmic function of the 
initial condition Eq. (|45|) . 

The generalization of this construction for higher order composite operators is possible, as well. The only condition 
is that the interaction Lagrangian for the elementary fields must be expressible in terms of the composite fields. 
The difficult question, the issue of the choice of the interactions in the effective theory can be hidden in the initial 
conditions for the effective action for the composite fields. This latter can partially be taken over from the elementary 
theory without interactions (aspects of the statistics) and partially be generated by the evolution equation (aspects 
of the interactions) which provides an iterative improvement of the effective theory (cf. Section [ill C|> . 

V. SUMMARY 

The functional rcnormalization group strategy was studied in this paper in the internal space in the framework of the 
scalar 4 model where either the mass or the coupling strength can be evolved. The mass-controlled suppression of the 
amplitude of the fluctuations realizes the functional generalization of the Callan-Symanzik scheme for the bare Green 
functions. The appropriate split of the evolution equation into the sum of two pieces provides a non-perturbative 
BPHZ-like renormalization of the 1PI Green functions. 

The control of the amplitude of the fluctuations by the interactions leads to the evolution equation generated by 
the change of the coupling strength. This scheme was exploited in order to obtain the evolution equation for the 2PI 
effective action. 

Finally, the similarity of the evolution equations in these schemes was used to construct an effective bilocal theory 
for the connected propagator. 

The advantage of the method, compared to the traditional external space renormalization group schemes, lies in 
its flexibility in choosing the parameter to evolve. This feature allows one to trace the dependence on certain scales 
in a model possessing several scales or on any of the coupling constants. This latter is important when symmetries 
prevent us to use momentum-space cutoff. 
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APPENDIX A: NOTATION, CONVENTIONS 

The space-time and Fourier-space integrals are denned as 

and the Fourier transform of a translation invariant propagator G x . y is defined as 8p+ q $G p — G p . q . The composite 
index X = (x, y), proves to be useful, e.g. J x y A x>y B Xy y — J x AxBx- The integration is frequently shown as scalar 
product, f ■ g — J x fx9x> A : B = J x F x ^ y G x ,y- The repeated indices are automatically summed/integrated over 
unless written explicitly otherwise. 

The functional derivative is defined in d-dimensional lattice as 



8 1 d 

where <\> L is the lattice field variable where the factor a~ d is needed in order to satisfy the relation 



(A2) 



4r I f {<!>,) = /'&,)• (A3) 



vyy J x 

In an analogous manner we have 



k~ r w <A4) 



Our convention for the derivative with respect to a symmetrical matrix Q s ^ m is 

2^ {eg™;-" yn "' + aes°r ,ym1 ') x ^ Vl 
i a 

a2d ogi y, y- 



5gf™ I -1 §— X -y 



where Q non s v m is a non-symmetrical matrix. It allows us the use of the chain rule 

8 f SB UV 8 



8A X y Ju,v ^A x ^y 8B U V 



(A5) 



(A6) 



APPENDIX B: ONE-LOOP EFFECTIVE ACTION 

The evolution equation for the 2PI effective action requires the knowledge of initial conditions imposed at vanishing 
or infinitesimal coupling strength. The loop-expansion for path integral Eq. I|33|) has two non-vanishing orders for 
g = 0, therefore these orders will be retained below. 

i. O {h°) 

The generating functionals are given by 

f(°)[S] = s[MJ}} - \UA ■ J • MA + ^Tr g l ■ j, (bi) 

in the tree approximation where the stationary point 4>o[J\ satisfies the equation 

™ =J-MJ), (B2) 

0<P \<f,=MJ\ 



15 



and the source J is related to the Green function Q via 

Gx,y = 4>0 x[J}4>0 y[J], 

i.e. 



'JO x 



(B3) 



(B4) 



Notice that <f>o[J] is non- vanishing and is given by this equation in the phase with broken symmetry only. The 
elimination of the source J yields 



f = s[MJ}} 



1 



.r..r "I 2 Qx.X ^| {Qx,x) 



(B5) 



In the symmetrical phase without tree-level condensate r tree [Cf] = holds. 



2. O (h) 

We write (f> = <po + 77, expand the exponent of the integrand in powers of r\ keeping the terms up to O (?y 2 ) and find 



w^[j] = -s[MJ}} + \<MA ■ J ■ MJ] - \tt& in 



Wc define the connected two-point function G by 

5WW[J] 

Gx,y = 4>0 x[J]4>0 y[J] + I^Gx,y = ~~P\~~r = ^0 x[J]4>Q y[A + & 



- J 



\MJ\ 



which gives 



G, 



v 2 ° x ,y 



PS[4>) 



S 2 S[<f>] 



\MJ] 



(B6) 



(B7) 



V,x 



\<t>o[J] 



and 



f«[£] - r tree [G] + ^Tr InG- 1 + ^Tr 



G ■ 



PS[<f>] 



(B8) 



(B9) 



up to a ^-independent constant, c.f. Ref. |l3j . In order to avoid the problem of non-convexity, we assume that the 
vacuum has no condensate, 4>a[J] — 0, i.e. G = ttG. 



APPENDIX C: INVERSE OF SYMMETRIC 4-POINT FUNCTIONS 



The quadratic part of the free effective action is 



*{x,y),(u,v) 



— ( G X G X ( e ip( - y ^ u ^ +iq ^ v ^ x " 1 + e i P(y~ v )+ i 9( u - x )\ 



(CI) 



for the translation invariant case. Therefore we seek the inversion of the operators of the form 



with symmetries Ar x y \r UtV ) = A( UV )( xy } = At x y \r VtU -\ which require A P:Q = A- q .- p = A q<p . It is easy to see that 
A^ 1 cannot be written in the form l|C2(l . But such operators can be inverted when they are acting on the space of 
symmetric matrices. In fact, the relation 



if, S +6 S ) — ~ f A A-iU^M^+e**-*™*-.)) 

2 \ u x,u u y,v \ u x,v u y,u) — ^ J ^PiQ p,q y ' J 



(C3) 
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yields I p . q = 1, A p £ = l/A p , q and 



for the free two-particle propagator. 



Sp.q 1 = -xQp 1 G q S Gp,q=2GpG q 



(C4) 



APPENDIX D: EQUATION OF MOTION 

The Schwinger-Dyson equation is obtained in this Appendix for the two- and four-point 1PI functions. The absence 
of the condensate is assumed, in order to simplify the expressions. 



1. Elementary field variable 



The path integral 



yields the Schwinger-Dyson equation 



which can be written as 



by means of the identity 



For 



Z[j] = e w ^ = / Z%] e - S M+^' 



ss[4>] 



ssy>] 



,wm = 



1 = 



i_ 5 | jWM 



SW[ 3 ] 



,W[j] 



Eq. (|D3(1 yields the Schwinger-Dyson equation 

, 9 







G 



Gq x ■ 



<s 2 r[0\ _1 s 3 w[j] 



y,y SjySjySjy 



-(Go-i)x. 



The first functional derivative of this equation is 



5 2 W[0] 

SjxSjy 



— Gq x ,y — — I Go 



9 f n S 2 W[Q] 6 2 W[0] g 



6 A W[Q] 

SjzSjy 5j z 5j z 6 J z "° x ' z Sj z 5j z Sj z Sj y 



-„ G 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



(D6) 



(D7) 



for j — 0. We shall express this relation by means of the 1PI vertex functions. By multiplying with G 1 and the 
inverse propagator we first arrive at 



<5 2 r[o] 



— 4- s 



g. S 2 W[0] g f S 4 W[0] 5 2 r[0] 



2 ' Sj x 8j x 6 ./, Sj x Sj x SjxSj 



(D8) 



Iz ^Vz^Vy 



In order to express the last term on the right hand side in terms of the effective action, one acts with the operator 



6 _ S<p S S 2 W[j] S 
Sj Sj Scj) SjSj 5(f) 



5 2 T 



(D9) 
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on the identity 



and finds 



1 _ 5 2 W[j] 5 2 T[ Y ] 
SjSj 54>S(f) 

6 3 W[j] 6 2 T[(t>} f S 2 W[j] S 3 T[(t)} S 2 W[j] 



IJ x 3 J{z} 



(D10) 



(Dll) 



where the shorthand notation Jr-. — J z ^ z has been introduced. Multiplication by the second functional derivative 
of W[j] in the index £3 gives 



S 3 W\ 



8 2 W[j\ S 2 W[j] S 2 W[j] <5 3 r[</>] 



$3xi8j X2 5j X3 J {z} Sj Xl 5j Zl 5j X2 Sj Z2 Sj X3 Sj Z3 84> Zl 84> Z2 5(j) Z3 
A further derivation leads to 

s 4 w[j] _ r r s 2 w[j] s 2 w[j] s 2 w[j] s 2 w[j] <s 4 r[</>] 



(D12) 



Sj xl Sj X2 Sj X3 Sj, 



Xi J{ Z } 



Sj xl Sj zl 5j X2 Sj Z2 6j X3 5j Z3 Sj Xi 8j 



24 vy Z2 K^ Z j vy Zi 



| S 2 W[j] S 2 W[j] S 2 W[j] <5 3 r[</>] 5 2 W[j] S 2 W[j] S 3 T[ Y } 
djxjj^ Sj x Jj Z5 Sj Zl dj Z6 8(j} Zi 8<t>z 5 Hz e $jxJj Z2 Sj X3 Sj Z3 5(j) Zl 6(j) Z2 S(j) Z3 

+ (xi < — > x 3 ) + (xi < — > Xi) 



(D13) 



which finally allows us to write Eq. (|D8|) as 



<5 2 r[0] 



x +9 _ 5 (s^my 1 

0x ' y 2 ' ? 'U>vJ,., 



<5 2 r[o] 



<5 4 r[o] 



S 2 T[0}\ fS 2 T[0] 



x, Z2 

-1 



x, Z3 vy- Zl °JYz 2 vYz 3 vry 

' x, Z2 V 6o,<i " ) x.. = ' ' 



The path integral 



leads to the relations 



S 2 W[0] 
<5 3 VF[0] 

<% ^1,2^2^3,4^2 ^5,6 



2. Product of two field variables 



Z[J] = e 



V[<p]e 



-S[ip]+^-J- v 



(D15) 



<W[0] _ S 2 W[0] 
8\J\,2 Sji5j 2 



G\ 2 — Gv 1, 



= if, 



(1,2), (3,4), 



S 3 Z[0] 



— 5v l,2<?v 3,4<?v 5,6, 



i^(l,2),(3,4)^v 5,6 - if(l,2),(5,6)£v 3,4 — (3,4),(5,6)£v 1,2 



where the two-particle propagator 



K 



( s 2 f[g v ] 
\ sgsg 



(D16) 



(D17) 
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has been introduced. (For the sake of simplicity the space-time coordinates x n are replaced by the corresponding 
indices n.) The analogue of Eq. (|D12|I is 



s 3 w[j] _ s 2 f[g}\ s 2 r[g}\ s 2 f[g}\ s 3 f[g] 



5\JxMJxMJx s \ sgsg / \ sgsg / \ sgsg / sg Zl sg z jg z , 

2 1 2 2 2 3 \ / Xi,Zi \ / x 2 ,z 2 \ / x 3 ,z 3 

Therefore, we have 

G\ 2 = g v 1,2, 

G l,2,3,4 = #(1,2), (3,4) ~ G v (1,2), (3, 4) 
G l, 2,3,4,5,6 - ~- K (l,2),Z 1 ^(3,4),Z 2 ^(5,6),Z ; 



(D18) 



sg Zl sg Z2 sg Za 

— [K{l,2),{3,5)G&$ + #(1,2), (4,5)^3,6 + #(1,2), (3,6^^5 + ^(1,2), (4,6^3,5 + G4 3 4 ^G^g + G§ 3 4 5G4 g 
+ ( 2l,3,4,6^ f 2,5 + G^^ gGJ g + G^ 3 5 ^G^ + G^gggG^ + Gi ;4 5 gGj 3 + Gj ,4,5,6Gi :3 ] (D19) 

with g v x,y the free two-particle propagator. The first functional derivative of the Schwinger-Dyson equation can be 
written as 

^x.y — ^^z.y^/v x.x g J {^-(x.x) ,(x,z) Gv {x,x),{x,z))Gv z.y (D20) 

Notice that the leading order perturbative result is recovered for K = g v only. When the two-particle sector is 
non-trivial by the choice of the independent variable g then the one-particle sector of the theory is naturally modified. 
In a similar manner the third functional derivative of the Schwinger-Dyson equation reads as 

9 f„ fr, r, rr ^v] 

■ Z 3 



K(x,y),(u,v) — £v (x,y),(u,v) + ^ J Gq x ,z (^K{z,z),Z l K{z,y),Z 2 K(u,v), 



' sg Zl sg Z2 sg Z} 

~^~K(z,z),(z,u)Gv y,v ~\~ K(z,z),(z,v)Gv y,u K(z,z),(u,v)Gv z,y K(z.y),(u,v)Gv z,z ~\~ z.zOv (z,y),(u,v) 

(D21) 



